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The two basic changes that were done to the quasi -one dimensional 
flow prog* am, where the replacement of the Runge-Kutta subroutine with a 
subroutine using a modified divided difference form of the Adams Pece 
Method and the replacement of the matrix inversion routine with a 
pseudo-inverse routine. The following system of first-order differ 
entlai equations Is solved In the program: 

A ^ « B, where A Is a n x n matrix, 


y and B are n-component column vectors. 


At the throat of the duct, If the flow Is approaching choked 
conditions then detA -*• 0 and the system Is singular. The pseudo-inverst^ 
routine determines, towards V ' solution of the linear system AZ » B, J 
the singular value decomposition, A = USV T . Householder biodiago- 
nal ization and a variant of the QR algorithm are used. The output from 
this routine Is S, V, and G, where G = U^B; S is diagonal and V,U are 


orthcg n nal; S also determines whether A Is singular. If one or more of 


the elements of S Is zero then A Is singular. In the subroutine, if 
S n < 10" 14 then S^" 1 is set equal to zero. The program then solves 


the folk i q system 
T dy _ ,.T, 


S V 




/ T & 
1 ax- 


il b 


S _1 G 


-1, 


Ml U 1 

(NAS A-CR-1L8691) NONLINEAR HAVE CHOKED 

INLETS Se*iann*id 1 Proqress Report, 1 Oct. 

1978 - 30 Apr. 1979 (Virginia Polytechnic 
Inst, and State Oniv.) 9 p HC A02/f*E rOI 

CSC L 20C 03/34 



Onclai 

22260 


* 


G 1 

If * — > 100 the program stops since this Indicates that y' Is becoming 
*11 

too large for the Integration routine to handle. Therefore Incompatible 
reflection coefficients (c R 's) have been chosen, If the above check 
Is passed, then y' Is Integrated. Solving for y' we have 

y’ « VS^G 

Now these values are passed to the Integration routine. Unlike the 
Runge-Kutta routine which would Integrate to the end of the duct even if 
Incorrect c R 's were Input, the present routine will not pass through the 
throat. The first cases considered were for one reflection coefficient 
where 

Pn(0) * c r Pm(0), and c R p?,(0) 

were specified. The initial conditions at x * 0 for the higher har- 
monics were set equal to zero. 

By varying c R the following plot was obtained showing the region of 
correct c R 's. This region was determined using two harmonics. Later 
cases using four and ten harmonics Indicated that the region did not 
shrink by much. 

The following examples Cjmpare the original Runge-Kutta version of 
the program with the modified version with the Adams-Pece and pseudo- 
inverse subroutines. The input to both runs is Identical except for the 
c R 's. For the Runge-Kutta case, c R = (0.,0.) and the results are for 
two cases, two and four harmonics. It is seen that as the throat is 
approached detA -*■ 0 and there is a substantial growth in the amplitudes 
of the harmonics. In the case of four harmonics there are very large 
jumps in the amplitudes near the throat. This indicates that incorrect 
c R 's were chosen. 





The example using the modified program was run with * (-.5, 

-3.3). As Is seen there are no large oscillations In the harmonics near 
the throat. This Implies that the Integration to the right end Is valid. 

I In the former cases, the Runge-Kutta routine could not detect the 

singularity and was Integrating across It. Therefore, the values obtained 
I by the Runge-Kutta method In the vicinity of the throat to the right end 

are Incorrect. If this occurs In the new version. It will terminate at 
the throat where 
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